Exact inhomogeneous Einstein-Maxwell-Dilaton cosmologies 
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We present solution generating techniques which permit to construct exact inhomogeneous and 
anisotropic cosmological solutions to a four-dimensional low energy limit of string theory containing 
non-minimally interacting electromagnetic and dilaton fields. Some explicit homogeneous and in- 
homogeneous cosmological solutions are constructed. For example, inhomogeneous exact solutions 
presenting Gowdy - type EMD universe are obtained. The asymptotic behaviour of the solutions is 
investigated. The asymptotic form of the metric near the initial singularity has a spatially varying 
Kasner form. The character of the space-time singularities is discussed. The late evolution of the 
solutions is described by a background homogeneous and anisotropic universe filled with weakly 
interacting gravitational, dilatonic and electromagnetic waves. 
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I. INTRODUCTION 

In the last decade string cosmologies attracted large 
amount of interest (see e.g. and references therein). In 
the traditional approach the present (low energy) string 
cosmology is in fact the classical cosmology where gen- 
eral relativity is generalized by including additional (in 
most cases) massless scalar fields. In general, one expects 
that the inclusion of these extra matter degrees of free- 
dom with the corresponding physical interpretation, may 
somehow resolve the long standing problems in cosmol- 
ogy. A recent interesting and promising development in 
this direction is the so called pre-big bang scenario 0. In 
the framework of this scenario one assumes that the ini- 
tial state of the universe is characterized by small string 
coupling and small curvature. This leads to an inflation- 
ary phase for sufficiently homogeneous initial conditions. 
There are, however, no natural reasons for the early uni- 
verse not be inhomogeneous. At present it is not clear 
how the large inhomogeneities may influence the pre-big 
bang scenario. In a more general setting, the question 
of whether our universe (which looks now homogeneous 
and isotropic) may arise out of generic initial data still 
lacks a complete answer in both general relativistic and 
string cosmology. 

That is why the construction and study of exact inho- 
mogeneous and anisotropic string-cosmological solutions 
are still of great importance. 

The inhomogeneous string cosmologies have been stud- 
ied by a number of authors. In [0], Barrow and Kunzc 
have presented inhomogeneous and anisotropic cosmo- 
logical solutions of a low energy string theory containing 
dilaton and axion fields when the space-time metric pos- 
sesses cylindrical symmetry . Their solutions describe 
ever-expanding universes with an initial curvature singu- 
larity. The asymptotic form of the solutions near the ini- 
tial singularity has a spatially varying Kasner-like form. 



In Fienstein, Lazkoz, and Vazquez - Mozo have pre- 
sented an algorithm for constructing exact solutions in 
string cosmology for heterotic and type - IIB superstrings 
in four dimensions. They have also presented and dis- 
cussed some properties of an inhomogeneous string cos- 
mology with S 3 topology of the spatial sections. Clancy 
et al. || have derived families of anisotropic and inhomo- 
geneous string cosmologies containing nontrivial dilaton 
and axion fields by applying the global symmetries of the 
string-effective action to a generalized Eistein- Rosen met- 
ric. Lazkoz H has presented an algorithm for generating 
families of inhomogeneous space-times with a massless 
scalar field. New solutions to Einstein - massless scalar 
field equations having single isometry, have been gener- 
ated in H by breaking the homogeneity of massless scalar 
field G2-models along one direction. 

It should be noted, however, that the inhomogeneous 
cosmologies in the framework of general relativity with 
certain physical fields as a matter source have been in- 
vestigated long before the time of the string cosmology. 
Exact inhomogeneous vacuum Einstein cosmologies with 
S 1 x S 1 x S 1 , S 2 x S 1 and S" 3 topology of spatial sec- 
tions have been found by Gowdy _M and studied after- 
wards by Berger || and Misner |0]7 Exact stiff perfect 
fluid inhomogeneous cosmologies have been studied by 
Wainwright, Ince and Marshma n JlO| ] . Later Charach 
fllT[ and then Charach and Malin |12[ (see also fDlfl) have 
found and studied exact inhomogeneous cosmological so- 
lutions to the Einstein equations with an electromagnetic 
and minimally coupled scalar field. These solutions have 
been interpreted as an inhomogeneous universe filled with 
gravitational, scalar and electromagnetic waves. 

The purpose of this paper is to present sufficiently gen- 
eral solution generating techniques which permit to con- 
struct exact inhomogeneous and anisotropic solutions to 
the equations of low energy string theory containing non- 
minimally coupled dilaton and electromagnetic fields - 



1 



the so called Einstein-Maxwell-dilaton (EMD) gravity. 
Examples of both homogeneous and inhomogeneous ex- 
act solutions will be also presented and their asymptotics 
will be investigated. 

Besides the fact that the EMD cosmologies are inter- 
esting in their own, there are at least two more moti- 
vations for considering the EMD cosmologies. First, a 
long standing question in cosmology is the existence of 
a primordial magnetic field. The string theory predic- 
tion of the non-minimal coupling of the dilaton to the 
Maxwell field gives an opportunity to study the problem 
in a more general dynamical context. Recently, in the 
homogeneous case this question has been addressed in 

B 

The second motivation is to investigate the character of 
initial cosmological singularities in the presence of matter 
fields arising from the low energy super-string theory. Ac- 
cording to the Belinskii-Khalatnikov-Lifhsitz (BKL) con- 
jecture the dynamics of nearby observers would decouple 
near the singularity for different spatial points. A vac- 
uum spatially homogeneous space-time singularity (IX 
Bianchi type) is described by BKL as an infinite sequence 
of Kasner epochs (" oscillatory or mixmaster behaviour" ) 
p5[ , p6| . BKL further speculated that a generic singu- 
larity should exhibit such a local oscillatory behaviour. 

Interesting special cases are the so-called asymp- 
totic velocity-term dominated (AVTD) singularities p7[ . 
Their characteristic feature is the fact that the spatial 
derivative terms in the field equations are negligible suf- 
ficiently close to the singularities. Such singularities are 
described by only one Kasner epoch, not by an infinite 
sequence of Kasner epochs i.e. an oscillatory behaviour 
does not exist. We think that it is important to study the 
behaviour of the singularities in the presence of matter 
fields coming form the super-string theory. In particu- 
lar, it is interesting to investigate the influence of non- 
minimal exponential coupling of the dilaton to the £7(1)- 
field on the character of the singularities. As it has been 
shown in Jl8| , the minimally coupled scalar field can 
suppress the oscillatory behaviour. 

Very recently the influence of the exponential dilaton- 
electromagnetic coupling on the character of initial sin- 
gularities has been studied by Narita, Torii and Maeda in 
p0| . These authors have considered T 3 Gowdy cosmolo- 
gies in Einstein-Maxwell-dilaton-axion system and have 
shown using the Fushian algorithm that space-times in 
general have asymptotic velocity-term dominated singu- 
larities. Their results mean that the exponential coupling 
of the dilaton to the Maxwell field does not change the 
nature of the singularity. It should be also noted that 
the exact solutions found in the present work lead to the 
same conclusion. 

Although, in the present paper we shall not consider 
the BKL conjecture in details, we believe that the present 
work may serve as a good ground for further studying 
of the above mentioned questions in the framework of 
Einstein - Maxwell-dilaton-axion gravity. 



II. CONSTRUCTING SOLUTIONS WITH ONE 
KILLING VECTOR 

We consider EMD-gravity described by the action 
A= — [ UR~2dip A*dip-2e- 2 fF h* F) 

16-7T J 

where * is Hodge dual with respect to the space-time 
metric g, R is the Ricci scalar curvature, tp is the dilaton 
field and F = dA is the Maxwell two form. 

First we will examine the case of space-time admitting 
one space-like, hyper-surface orthogonal Killing vector 
X. Physically this situation corresponds to a universe 
with homogeneity broken along two space directions. In 
the presence of a Killing vector the Bianchi identity dF = 
and the Maxwell equations d * e~ 2<p F = give rise to 
the following local potentials ^> e and 

rf* e = -i x F , dV m = e- 2v i x * F. 

Here ix is the interior product of an arbitrary form 
with respect to X. 

The hyper-surface orthogonal Killing vector X natu- 
rally determines a three dimensional space-time subman- 
ifold £ with metric §l) 

P = Mg — X ® X 

where Af = g(X, X) is the norm of the Killing vector. 

We note that throughout the text we denote the Killing 
fields and their naturally corresponding one-forms by the 
same symbols. 

The projection on S of the Einstein equations reads 

K = 2dtp ®dip+ -^dM <E)dJ\f + 
2N 2 

^ (e- 2v d^ e <g> d* e + e 2lfi d^ m <g> dtf m ) (1) 

where 1Z is Ricci tensor with respect to the projection 
metric P. The Maxwell equations take the form 

S (e-^M^d^e) = , d) (tpe/f-id&m) = 0. (2) 

Here d) = *d* is the co-derivative operator. 
The projection of the Einstein equations along to the 
Killing vector gives 

AW - g{dJV, dJV) = 
-2Ne- 2 * (e- 2 ^g(d^ ei d^ e ) + e 2 ^ 5 (d* m , d* m )) . (3) 

Finally, the vanishing of the twist T = | * dX of the 
Killing field leads to the constraint 

dt> e A d* m = 0. (4) 
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The equations ([!]) , (0) , @ and ([|) are equivalent to the 
EMD - gravity equations in the presence of a spice-like, 
hyper-surface orthogonal Killing vector. 

The constraint ([|) implies that the one-forms d^f m and 
dty e are proportional i.e. there exists a function T such 
that d^m = T d^ e . Here we will consider the simplest 
case in which one of the potentials \& vanishes, say \& e = 
0. 

In the case \l/ e = the full system dimensionally re- 
duced equations can be obtained from the following ac- 
tion 



(5) 



Now we introduce the following symmetric matrix S G 
GL+(2,R): 



N + 2^ 2 e 2 ^ v / 2* m e 2 ^ 



5' 



2<p 



Making use of the matrix 5*, we may write the action 
(||) in the form 



(6) 



The action (|J) has GL(2, R) as a group of global sym- 
metry for fixed projection metric P. Explicitly the group 
GL(2, R) acts as follows : 

S -> ASA T 
where A £ GL(2,R). 

Note, that in the case when one of the electromag- 
netic potentials vanishes the dimensionally reduced EMD 
equations can be viewed as three-dimensional Einstein 
gravity coupled to a nonlinear a - model on the factor 
GL(2,R)/0(2,R), 

The GL{2 1 R) symmetry may be employed for gener- 
ating new exact solutions from known ones. In particu- 
lar, it will be very useful to employ the GL(2, R) sym- 
metry to generate exact solutions with nontrivial elec- 
tromagnetic field from any given solution to the vac- 
uum field equations (i.e. pure Einstein equations) or 
the dilaton-vacuum equations (Einstein equations plus 
a dilaton field). 

For example, starting with a dilaton-vacuum seed so- 
lution 



the nonlinear action of GL(2, R) gives 



N = (detA) 



2 Kd e 2y - d 
c 2 N v d + d 2 e 2 ^« d 



3 2{p = c 2 Af vd + d 2 e 2ip " d , 



*m = 



1 acN vd + bde 2v - d 
71 c 2 N vd + d 2 e 2 ^ 



where 



A = 



a b 
c d 



It is worth noting, however, that some of the GL(2, R)- 
transformations are pure electromagnetic gauge or rescal- 
ing of the solutions. Only the 0(2, ^-transformations 
lead to essentially new solutions. 

Exact solutions to the EMD gravity equations can be 
constructed from vacuum solutions by the method de- 
scribed in (22), too. In order to obtain solutions by 
this method one assumes that the matrix S depends on 
space-time coordinates through a harmonic potential f2 
(D % Di fl — 0). Requiring then the constraint 



_i fdS ds- 1 

4 r \dfl dtt 



= 1 



to be satisfied, the EMD equations are reduced to the 
vacuum Einstein equations 



Hi 



2DinDjQ, , 
D*Din = 



and a separated matrix equation 



A. f q- 1 — 
dn [ dQ 



Therefore, for every solution to the vacuum Einstein 
equations, the solutions of the matrix equation give 
classes of exact EMD gravity solutions (see for details 

H)- 

It should be noted that the seed solutions used for the 
construction of exact EMD solutions via the above de- 
scribed methods, must admit at least one Killing vector. 
The case with two commuting hyper-surface orthogonal 
Killing vectors will be considered in the following section. 



K,d 
e 2 <^ 
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III. CONSTRUCTING SOLUTIONS WITH TWO 
KILLING VECTORS 

A. General equations 

In this section we consider the EMD gravity equations 
when the space-time admits two commuting space-like, 
hyper-surface orthogonal Killing vectors X and Y. In 
this case the metric can be written in the form 



ds 2 = e 2h-2 P 



(dz 2 -d£ 2 ) + e 2p dx 2 + g e~ 2p dy 2 (7) 



where h, p and g are unknown functions of z and £ 
only. Therefore the Killing vectors are given by X = 



Y = 



' hi ' 



The physical properties of the metric (0) depend on the 
gradient of the norm g of the two-form X f\Y . The case 
corresponding to d^g being globally space- like or null de- 
scribes cylindrical and plane gravitational waves. Metrics 
where the sign of d^gd^g can change, describe colliding 
plane waves or cosmological models with space-like and 
time-like singularities. Metrics with globally time-like 
d^g describe cosmological models with space-like singu- 
larities. 

In the present paper we shall consider only the globally 
time-like case g — £ 2 . It should be noted that the choice 
g = £ 2 is dynamically consistent with the EMD gravity 
equations because the Ricci tensor R satisfies 

g(Y, Y)R(X, X) + g(X, X)R(Y, Y) = 0. 

In the presence of a second Killing vector the dimen- 
sional reduction can be further continued. Without going 
into details we directly present the dimensionally reduced 
equations 

%P + - 8 2 p = e- 2 *- 2 ? ((9 z * e )2 _ (d^ e ) 2 ) 

+ e 2 ^ 2p ((d z ^ m ) 2 -(d^ m ) 2 ), (8) 



%<p + -%> - %<p = e- 2 *- 2 ? ((d z ^ e ) 2 - (d^ e ) 2 ) 

- e 2 «>- 2p ((0,* TO ) a - (d^ m f) , (9) 

d^ed z *m = d z V e d^ m , (10) 
B e (^- 2 Pd^ e ) - d z (ee-^- 2p a z * e ) = o, (11) 
flfe (Ce 2 ^- 2 Pd^ m ) - d z {Ze 2 v- 2 *>d z * m ) = 0, (12) 

= [d^f + (d z p) 2 + {d iP f + (d zP ) 2 + 

e -2 V -2 P ( (9 ^ e )2 + ( d ^ e ) 2 ) + 
e -2 V -2 P ( (9 ^ m) 2 + ( 5 ^ m )2) , (13) 



-d z h = 2d^(pd z (p + 2d(j>d z p ■ 
2e- 2 *- 2p d^ e d z * e + 2e 2 *- 2p d^ m d z *„ 



(14) 



Here we should make some comments. We have written 
the system of partial differential equations in terms of the 
electromagnetic potentials fy e and \t m . These potentials 
are derived with respect to the Killing vector X. In the 
same way we may use the electromagnetic potentials with 
respect to the Killing vector Y denoted respectively as 
^ and Moreover, we may also use a mixed pair 

of potentials, say ^ e and W^f . In general, the transition 
between different pairs of potentials may be performed 
by the following formulas 



-2ip-2p 



d z V e 



d z ^ Y e =-ie 2 ^~ 2p d^ r 



£e 2 *- 2p d z * r 



For the reader's convenience we have presented in the 
Appendix the system (||) - (14) written in terms of the 
mixed pair f e = uj and = \ ■ 

Let us go back again to the system of nonlinear partial 
differential equations (|J) - (|T4|). We have a complicated 
system of coupled nonlinear partial differential equations 
and it seems that finding all its solutions is a hopeless 
task. Nevertheless, as we will see a large enough class of 
solutions can be found. 



B. Solution generating method 

Here we consider the case in which one of the electro- 
magnetic potentials vanishes. For definiteness we take 
$ e = 0. Let us introduce the new potentials u — p — ip, 
4> = p + ip, = V2~y m and h s = 2h. Then the system 
(^) - ( |l4| ) may be rewritten in terms of the new variables 
as follows 



%4> + jdtf - d 2 <f> = o, 

din + iflfet, - d 2 u = er 2 - {(d z ^ m ) 2 - (d^ s m ) 2 ) , 
«£ (tie- 2u d^ s m ) - 3 Z (£e- 2u d z y s m ) = 0, (15) 



^h s = (d^) 2 + {d^f + (dtp 3 ) 2 + (d zP s ) 2 + 

e~ 2 « ((d^ s m ) 2 + (d z * s m ) 2 ) , 
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<P = (f> +/3 HO, 



(18) 



O z h s = 28^8^ + 2d iU d z u + 2e- 2u d^f n d z ^! i r 



It is not difficult to recognize that the system (|l5| ) 
coincides with the corresponding one for the Einstein- 
Maxwell gravity with a minimally coupled scalar field <j) 
for the metric 



ds 2 



„2fc°-2ii 



{dz 2 - di 2 ) + 



z 2u dx 2 +e~ 2u edy 2 . (16) 



Here is the corresponding non-vanishing electro- 
magnetic potential in the case of a minimally coupled 
scalar field. 

In this way we have proven the following proposition 

Proposition B.l Let u, <j>, h s and be a solution to 
the Einstein- Maxwell equations with a minimally coupled 
scalar field for the metric (|7(jj. Then p = ^(u + <fi). 
ip = \{4>— u), h = T;h s and ty m = -^&m form a solution 
to the equations of EMD gravity for the metric Q). 

This result allows us to generate exact solutions to the 
EMD gravity equations for the metric (|t|) directly from 
known solutions of the Einstein - Maxwell equations with 
a minimally coupled scalar field. 

Before going further we shall make some comments. 
In the case when the space-time admits two space-like 
Killing symmetries and one of the electromagnetic po- 
tentials is taken to vanish, it may be shown that the 
equations governing the transverse part of the gravita- 
tional field, dilaton and non-vanishing electromagnetic 
potential are written in the compact matrix form 



(17) 



The equation ( [l7| ) is just the chiral equation. There 
are powerful solitonic techniques for solving this equation 
(see p3[). Although, the study of cosmological soliton 
solutions in EMD gravity seems to be very interesting, 
we will not consider them in the preset work. 



In 



*fr = l tanh( ^ + ^ m(e))) 



where (f>o, </>i, a$ and /3 are arbitrary constants. 



(19) 



(20) 



The above seed solution (18) - ( pOj ) gives the follow- 
ing homogeneous solution to the EMD equations. The 
potential tyj is 



*5 = 



l 



31 , a 



The metric has the form 

ds 2 = A 2 ((,)(dz 2 - d£ 2 ) + B 2 ((,)dx 2 + C 2 {£)dy 2 
where 



D 



> + e 2 £ 2 



C 



!(0 = e -^o£i-/3o (e^^+e^C^) 



Here 70 is a constant. 

The asymptotic behaviour of the expansion factors and 
the dilaton field in the limit £ — > is as follows 



IV. EXAMPLES OF EXACT SOLUTIONS 

A. Homogeneous solutions 

Homogeneous solutions are obtained when we have de- 
pendence only on the "time - coordinate" £. In this case 
the differential constraint reduces to 8^ e d^ m = 0. 
Thus we may choose one of the potentials to vanish. 
Therefore we may apply the proposition from the pre- 
vious section to generate homogeneous solutions starting 
from the corresponding homogeneous solutions for the 
case of a minimally coupled scalar field. In our notations 
the homogeneous solutions for a minimally coupled scalar 
field are l|13f 



c 2 (o~e- po+l -^, 

V~~(l-A)-^|ao |) MO- 

The asymptotic form of the expansion factors and the 
dilaton field in the limit £ ^> 1 is 
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B. Inhomogeneous solutions 



3 2 (£)~£ 1+/3o + ^. 



V~-~(l-A> + ~ \<*o I) MO- 

Introducing the scynhroneous time rfr = the 
line element in the limits £ — > and £ 3> 1 can be written 
in the Kasner form 



Here the proposition (B.l) will be applied again, this 
time to generate inhomogeneous exact solutions to the 
EMD equations. As a seed family of solutions we take the 
Charach family of solutions describing inhomogeneous 
cosmologies with minimally coupled scalar and electro- 
magnetic fields and with S 1 x S 1 x S 1 topology of the 
spatial sections. The Charach family of inhomogeneous 
solutions and their asymptotics in the two limiting cases 
are presented in the Appendix. 



According to the proposition (B.l) the inhomogeneous 
EMD metric is given by 



ds 2 



-dr 2 + r^dx 2 + T 2 ^-dy 2 + r 2 ^dz 2 



ds l 



where 



^ h ~P\dz 2 - dt 2 ) + e 2 Pdx 2 + ee~ 2p dy 2 



where the Kasner exponents are defined by 



Pi 



Pi 



1 + 0o T \ I «o 



i(lT|ao|) 2 + (i-/3o) 2 + | ' 



l - A) ± h 



i(l T |aol) 2 + (i-/3o) 2 + | ' 



2p=ln(2£cosh(~<£)) + </>, 

h = ~ln(£)+rn(2cosh(^)) + 
-F(<j> ,a , A n , B n ; z) + F(<f> ,f3 , C n , D n ; £, z). 
The dilaton field is 



P3 



KIT I a |) : 



■A) 



,2_ 1 
' 2 



i(l T |ao|) 2 + (i-/3o) 2 + 



3 ' 



as the sing " — " refers to the limit £ — > while the sign 
f " is for the limit £ > 1. 
The dilaton filed is given by 



tp = -(f>- - In ( 2£cosh(-0) 



The inhomogeneous cosmological solutions introduce 
characteristic length scales. In fact each normal mode 
has its own characteristic scale. The horizon distance in 



the "z" direction is given by ds 2 



and hence 



where 



(p ~ a ln(r) 



1 - A) T \ot 



i(l T |ao|) 2 + (i-/3o) 2 + 



The parameters px, pi and p% satisfy the Belinskii 
Khalatnikov relations 

Pi +P2 +P3 = 1, 



pi 



-pi 



-pl = l-2a 2 



Sz= d£_=^. 



In this way n£ can be viewed as the ratio of the hori- 
zon distance in the " z" direction to the coordinate wave- 
length A„ i.e. n£_ — j^. 

There are two limiting cases to be considered. The 
first case is when the wavelength is much larger than the 
horizon scale (n£ <C 1). The second case is when the 
wavelength is much less than the horizon scale (n£ 3> 1). 

In the first case (n£ < 1) the asymptotic form of the 
EMD metric is 



ds 2 = A 2 (£, z)(dz 2 - d£ 2 ) + B 2 (i, z)dx 2 + C 2 (£, z)dy 2 
where 
A 2 (t,z) = 
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IS 



e -4>,{z)g.-l3(z) ^Mz)^a(z) + e -^.W^-5«Wj _1 . 

In the limit n(«l the asymptotic form of the dilaton 

^~ \ [4>*{z) + sign(a(z))4>*(z)j - 

\{l-P{z)-\\<*{z) 

In order to discuss the cosmological solutions near the 
singularity we have to consider homogeneous limit for £ 
approaching zero. Introducing the proper time, in the 
homogeneous limit, by 

r = J A(Odt 

we find that the metric has the Kasner form 

fl^~(-l,r 2 " 1 ,r 2 »,r 2 '"). 
Here the Kasner indexes are spatially varying 



Pi(z) = T 



Vi(z) = - 



Ps(z) 





l + /3(z) + || 


a(z)\ 




|a(*) + (I 






l-/?(z)±± | 


a(z) | 




\a(z) + (I 


-mr + i 




1 I) 2 + (i 


-P{z)f-\ 




1 I) 2 + {\ 


-P(z)) 2 + I 



and satisfy the Belinski- Khalatnikov relations 

Pi(z) +p 2 (z) + p 3 (+) = 1, 
p 2 1 (z)+p 2 2 (z)+pl(z) = l-2a 2 (z), 



where 



l-/3(z)Th<*(z) 



2,3- 



\{l T \a{z)\y + {\-P{z)y + 

In the second case when n£ 3> 1 the asymptotic form 
of the EMD metric is as follows. 

The first sub-case is for ao = 0. Then the EMD metric 
is given by 



where 



and 



h = diag (o, l=e +f3o H(t, z),-^=e- 0o m, z), 
When a ^ we have 

r\ = diag (-£Ki+«o) 2 -K/3o-|)M e i^ £i+tH«o+/3o 5 

t l-i|a |-/3o ti(l+ao) 2 + (/3o-^ 2 - 1 1 



V+a )-+(Po--i)---j e iKZ 



h = diag(o, (gfc z) + ±sign(a )H(Z, z)) 

- £^^ Z ) + I si5 n(a )i?(C, *)), o) . 



The asymptotic form of the dilaton in the limit n£ 3> 1 , 
for both «o = and «o 7^ 0, is given by 



I ln(0 + 



z) + isz ff n(a )iJ(£, z) 



The explicit form of the solutions and their asymp- 
totics allow us to make some conclusions. The obtained 
solutions fall in the category of asymptotic velocity-term 
dominated space-times. When the singularity is ap- 
proached the spatial derivatives become negligible with 
comparison to the time derivatives. Sufficiently close to 
the singularity the evolution at different spatial points 
is decoupled and the metric is locally Kasner with spa- 
tially dependent Kasner indices satisfying the Belinskii- 
Khalatnikov relations. 

The late evolution of the exact cosmological solu- 
tions is described by a homogeneous, anisotropic universe 
with gravitational, scalar (dilatonic) and electromag- 
netic waves. The non-minimal dilaton-electromagnetic 
exponential coupling influences mainly the homogeneous 
background universe rather the scalar and electromag- 
netic waves on that background. The former may be 
considered as minimally coupled scalar and electromag- 
netic waves up to higher orders of |. 
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V. CONCLUSION 
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In this paper wc have shown that it is possible to 
find exact inhomogeneous and anisotropic cosmologi- 
cal solutions of low energy string theory containing 
non-minimally interacting dilaton and Maxwell fields - 
Einstein-Maxwell -dilaton gravity. First we have consid- 
ered space-times admitting one hyper-surface orthogonal 
Killing vector. It has been shown that in the case of 
one vanishing electromagnetic potential the dimension- 
ally reduced equations possess a group of global symme- 
try GL(2,R). We have described an algorithm for gen- 
erating exact cosmological EMD solutions starting from 
vacuum and dilaton-vacuum backgrounds by employing 
the nonlinear action of the global symmetry group. This 
algorithm is especially useful for generating exact cos- 
mological EMD solutions with only one Killing vector, 
starting with Gi-dilaton-vacuum background. Another 
method which permits to construct exact EMD solutions 
starting from solutions of the pure Einstein equations has 
been briefly discussed, too. 

In the case when the space-time admits two commut- 
ing, hyper-surface orthogonal Killing vectors we have 
given a method which allows exact inhomogeneous cos- 
mological solutions to the EMD equations to be gener- 
ated from the corresponding solutions of the Einstein- 
Maxwell equations with a minimally coupled scalar field. 
Using this method, as a particular case, we have obtained 
exact cosmological homogeneous and inhomogeneous so- 
lutions to the EMD equations. The initial evolution de- 
scribed by these solutions is of a spatially varying Kasner 
form. The intermediate stage of evolution occurs when 
the characteristic scales of the inhomogeneities approach 
the scale of the particle horizon. This stage is character- 
ized by strongly interacting non-linear gravitational, dila- 
tonic and electromagnetic waves. The late evolution of 
the cosmological solutions is described by a background 
homogeneous and anisotropic universe filled with weakly 
interacting gravitational, dilatonic and electromagnetic 
waves. 

The cosmological solutions found in the present work 
fall in the category of AVTD space-times. Near the sin- 
gularity the dynamics at different spatial points decou- 
ples and the metric has a spatially varying Kasner form. 
Therefore, the non-minimal coupling of the dilaton to the 
Maxwell field does not change the nature of the singular- 
ity. The solutions are sufficiently generic and therefore 
we should expect that the T 3 space-times in EMD gravity 
have AVTD singularities in general. 

Finally we believe that the present paper will be a suffi- 
ciently good background for studying similar problems in 
the more-general case of Enstein-Maxwell-dilaton-axion 
gravity. 
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APPENDIX A: REDUCED EMD SYSTEM IN 
TERMS OF THE POTENTIALS * B = cu AND 

%p + \d iP - dip = e- 2 *- 2 ? {{d z u) 2 - {d^f) 
+ 1 e- 2 ^ ((da) 2 (d zX f) 

8 2 V + ±dz<p - 2 z cp = e- 2 ^ {{d z uj) 2 - (d^) 2 ) 
^ 2 e- 2 ^((d a ) 2 -(d zX ) 2 ) 

d^ujd^x = d z ujd z x 

d 2 oj + — d z u> = 
2(d^ip + d^p)d^u> - 2(d z ip + d z p)d z ui 

dfx ~ 7%X - d 2 z X = 
2(d^ - d^dt-x - 2(<9 Z ^ - d z p)d z x 

jdth = (d^f + (d zV f + (dtp) 2 + (d zP ) 2 + 
±e- 2 ^ ((d,xf + (d 5X ) 2 ) 

^d z h = 2d^fd z (p + 2d^pd z p + 
2 e -*p-*Pd i u>d z w + 2e- 2v+2p d a d zX 
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APPENDIX B: CHARACH FAMILY OF EXACT 
SOLUTIONS 

For the reader's convenience we present here the 
Charach family of solutions written in our notations. 
The minimally coupled scalar filed is given by 

(p=(j> o +0o HO + 

OO 

(CnJoH) + D n N (n£)) cos(n(z - z n )). 

n=l 

Here 4>o, Po, C n , D n are arbitrary constants and Jq(.) 
and Nq(.) are respectively the Bessel and Neumann func- 
tions. 

The metric function u and the non-vanishing electro- 
magnetic potential Y are correspondingly 



a = In ( 2£cosh(i</>) ) . 



*f Y = const + i tanh(^) 

where <fi is an auxiliary scalar field which is of the form 

<f> = <f>o + a Q ln(£) + 

OO 

cos(n(z - Zn)). 

n=l 

The longitudinal part of the gravitational field is 
h s = ln(£) + 2 ln(2 cosh(i</>)) + 
-F(cf> , a ,A n ,B n ;£, z) + 2F(0 O , Po, C n , D n ; £, z) 
where F is a solution to the system 

^F=i((9^) 2 + (9 Z 0) 2 ), 



and 



d z F = d^4>d z 



For more details see |13fl . 

In the limit £ — * the minimally coupled scalar field 
and the auxiliary field behave as 

0~^,(z)+/?(z)ln(O, 

where 

2 °° 

a(z) = a + - B n cos(n(z - z n )), 

n=l 
2 °° 

P(z) = fa + - y~) D n cos(n(z - z n )) 



4>*(z) = 4>o + 



J2 (Cn + -A l (7 + ln(^))) cos(n(z-z n )), 



71=1 



4>*{z) = 4>o + 

( A n + - B n{l + ln(-)) J cos(n(z - z n )). 
Here 7 is the Euler constant. 

The asymptotic form of the Charach metric can be 
written as 



ds 2 = Alii, z){dz 2 - d( 2 ) + Bs(£, z)dx 2 + C 2 {t z)dy 2 
where 



B 



c s % z) = ( e 4*' w c* aW + e -5^w^-s«w^ 1 . 

The asymptotic behaviour of the scalar fields in the 
Charach solution in the high frequency regime (n£ S> 1) 
is 

0- 0o + ao lnC + r^-H"(C,2), 

where 



H((,z)=Re V (A w -tB N ) 



g-*(«n-f ) 

y/2n\n\' 



00 -i{ z _») 

H{Z,z)=Re J2 [C\n\-iD\n\)- -^N*+~) 



n— — oo,^0 
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The functions H(£,z), H(£,z) satisfy the D'Alembert 
equations 

d 2 H(tz)-d 2 z H(t,z)=Q, 
d 2 H(^z)-d 2 H(tz)=0. 

When ao = the asymptotic form of the metric is 
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where 



and 



h = diag \0, A^H (£, z)/4£, -^(6 z)/4£ 
The constant K is given by 
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